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Liquid state nuclear magnetic resonance (NMR) techniques have produced some spectacular suc- 
cesses in the construction of small quantum computers, and NMR is currently by far the leading 
technology for quantum computation. There are, however, a number of significant problems with 
any attempt to scale up the technology to produce computers of any useful size. While it is proba- 
ble that some of these will be successfully sidestepped during the next few years, it is unlikely that 
they will all be solved; thus current liquid state NMR techniques are unlikely to provide a viable 
technology for practical quantum computation. 



INTRODUCTION 

Nuclear Magnetic Resonance (NMR) @, §, |, |, | § 
is almost unique among potential quantum technologies 
in that it has already been used to build small quan- 
tum computers @, | |, 0, 0. Although other tech- 
niques have been used to implement quantum logic gates, 
such as the ion trap controlled-NOT gate NMR 
provided the first complete implementation of a quan- 
tum algorithm ]l0| . Since then progress has been ex- 
tremely rapid, with demonstrations of Dcutsch's algo- 
rithm |h], [r|, Grover's quantum search pi, ver- 
sions of the Deutsch- Jozsa algorithm with three [ |17| ] and 
five |l8| qubits, and quantum counting ]l9| . 

In addition to being one of the most successful 
quantum computing technologies, liquid state NMR is 
also among the oldest. Although explicit experimental 
demonstrations of NMR quantum computation only date 
from 1996 §, many "conventional" NMR experiments, 
such as COSY g [H) and INEPT || can in retro- 
spect be reinterpreted as quantum computations. These 
and related experiments, based on coherence transfer se- 
quences, have been in use all over the world for decades, 
and are regularly used to study complex biomoleculcs 
containing thousands of nuclear spins p3]. Even older 
than these are simpler experiments such as selective pop- 
ulation transfer pif , which corresponds to the implemen- 
tation of a controlled-NOT gate, although in this case the 
gate was not applied to spins in superposition states. 

The rapid progress made in NMR quantum compu- 
tation builds on this pre-existing experimental sophisti- 
cation: decades of experience in manipulating nuclear 
spins in coherent states has inevitably resulted in a wide 
range of "tricks of the trade" . Given this success it may 
seem surprising that some early papers suggested that 
this approach might be limited to quantum computers 
containing about 10 qubits |^5|, ^6|, while many more 
modern estimates are similar at 10-20 qubits. The expla- 
nation for this is simple: the highly developed nature of 
NMR experiments, in comparison with many other puta- 
tive quantum technologies, means that the limits of the 
technique are fairly well known and understood. Thus 



unduly optimistic predictions about the power of NMR 
quantum computation may be easily debunked by expe- 
rienced NMR spectroscopists. 

STRUCTURE AND SCOPE 

It is conventional when assessing proposals for quan- 
tum computation technologies to consider how the tech- 
nique can be used to implement five basic elements re- 
quired to build a quantum computer, and then to discuss 
whether these techniques can be scaled up for use in com- 
puters with large numbers of qubits. This approach is not 
appropriate for NMR, as it is already clear that all the 
five basic elements can be constructed; indeed they have 
all been experimentally demonstrated. Thus, the only 
topic remaining to be addressed is that of scaling. It is, 
however, useful to consider each of the five basic elements 
in turn, and discuss their impact upon the practicality of 
building a large NMR quantum computer. 

The sections below are largely concerned with what 
might be called "conventional" liquid state NMR quan- 
tum computers, by which I mean computers implemented 
using standard techniques from liquid state NMR with 
spin-i nuclei. I will, however, briefly discuss some rela- 
tively simple techniques, such as the use of optical pump- 
ing or the use of liquid crystal solvents, which may allow 
small extensions in the range accessible to NMR quan- 
tum computers without greatly altering the underlying 
physics. I will not consider solid-state NMR: although ul- 
timately based on the same physical interactions the solid 
state NMR Hamiltonian |27], |29| is much more com- 
plicated than the liquid state form, and liquid and solid 
state NMR form two largely separate sub-fields. While 
some solid state NMR experiments could be considered as 
implementations of quantum simulations (e.g., ]3(|), they 
have not as yet been used to build general purpose quan- 
tum computers. Similarly, while nuclei with spin quan- 
tum numbers greater than \ (quadrupolar nuclei) are of 
some theoretical interest |n], |3^], their use raises con- 
siderable experimental difficulties. Finally I will largely 
ignore proposed systems which use single isolated atomic 
nuclei in solid state devices, such as the proposal due 
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to Kane |53| [34j; although such systems are ultimately 
based upon NMR, they differ from current liquid state 
NMR implementations so profoundly that it is difficult 
to draw detailed parallels. 



INITIALIZATION 

Initialization is the process of placing a quantum com- 
puter in some well defined initial state, typically |0) = 
|000. • .0), prior to beginning the computation. In prin- 
ciple any initial state is as good as any other, but in prac- 
tice |0) is the most widely chosen, both because it corre- 
sponds to the traditional starting point of many quantum 
algorithms and because it often corresponds to the sys- 
tem's energetic ground state. This is indeed the case in 
NMR quantum computation, where the computational 
basis corresponds with the natural experimental basis. 

As initialization requires that the quantum computer 
be placed in the state |0), independent of its state before 
the beginning of the initialization process, it is clear that 
it cannot be achieved by any unitary process; thus ini- 
tialization schemes must be quite different from quantum 
logic gates. As the desired initial state is usually an ener- 
getic ground state, initialization is typically achieved by 
cooling. This is not a practical approach within NMR, 
as the energy gaps involved are tiny compared with the 
Boltzmann factor at room temperature. 

The energy gap between nuclear spin levels in NMR ex- 
periments is principally determined by the Zecman inter- 
action between the nucleus and the applied magnetic field 
(with the exception of the quadrupolar interaction, which 
does not occur in spin-i nuclei, all other nuclear spin 
interactions are small compared with the main Zeeman 
interaction). The Zeeman splitting, AE = hv = fvyB, 
is usually described in terms of the corresponding Lar- 
mor frequency, z/, and is proportional to the magnetic 
field strength, B, and the gyromagnetic ratio, 7, which 
is an intrinsic property of the nucleus. For the magnetic 
field strengths typically used in NMR experiments (2.3- 
21.1 T), the Larmor frequency of 1 H nuclei lies in the 
range 100-900 MHz, corresponding to an energy of 0.4- 
3.7 /xeV. For all other nuclei (with the exception of the 
radioactive nucleus tritium, 3 H), 7 is lower than for 1 H, 
with a corresponding reduction in AE. 

These energies are much smaller than the Boltzmann 
energy at room temperature (about 25meV), and so at 
thermal equilibrium the excess population in the lower 
Zeeman level is tiny, less than one part in 10 4 . For this 
reason conventional liquid state NMR was long ruled out 
as a practical technology for quantum computation. In 
1996, however, it was realised that it is not strictly nec- 
essary to start quantum computations from a pure state: 
a pseudo-pure, or effective pure state will suffice R]. 



Pseudo-pure states 

A pseudo-pure state is a mixed quantum state corre- 
sponding to a mixture of the desired quantum state, 
and the maximally mixed state, 1/N, where N = 2 n is 
the dimension of the Hilbert space describing the sys- 
tem. To perform a quantum computation it suffices to 
form the pseudo-pure ground state, 

p=(l-e)l/iV + e |0)(0|. (1) 

An otherwise error-free quantum computer which begins 
a computation in such a state will end up in the state 

p' = (l-e)l/JV + e|V)<V| (2) 

where is the result of the computation (assumed for 
the moment to be an eigenstate in the natural basis). 
This result may be immediately deduced by noting that 
any quantum computation corresponds to some unitary 
evolution of the quantum state; such evolutions are linear 
and have no effect upon 1. 

A quantum computer in this mixed state will return 
the correct answer with probability e + (1 — e)/N, and a 
wrong answer with probability (N — 1)(1 — e)/N (note 
that the maximally mixed state itself contains a frac- 
tion 1/N corresponding to the ground state). Clearly it 
is possible to determine the desired answer by statisti- 
cal analysis of the results of a sufficiently large number 
of repetitions of the computation, where the number of 
repetitions required depends on e. Equivalently it is pos- 
sible to use an ensemble of quantum computers, and de- 
termine the ensemble averaged result of the computation; 
as long as the ensemble is sufficiently large this process 
will unambiguously point to the desired answer. 

This latter approach is precisely that adopted for NMR 
quantum computation, and indeed for conventional NMR 
spectroscopy. NMR transition frequencies are so low that 
it is essentially impossible to detect a single transition, 
and so it is necessary to use macroscopic samples, typ- 
ically containing about 10 17 molecules, with an excess 
population of about 10 13 nuclei in the low energy state. 
It might seem that the small signal from the excess popu- 
lation would be swamped by a huge background, but this 
is not the case as the maximally mixed state gives rise 
to no overall signal. This is because the signals from dif- 
ferent components in the maximally mixed state cancel 
each other out; the operators corresponding to the effec- 
tive observables in NMR spectra are all traceless. Thus 
the NMR signal arises entirely from the small excess pop- 
ulation, and the signal from an NMR quantum computer 
in the pseudo pure state, equation (Q), is identical to that 
from one in the pure state except that the signal 

strength is reduced by a factor of e. 
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Assembling pseudo-pure states 

While pseudo-pure states offer the theoretical possi- 
bility of performing quantum computations with mixed 
states, this approach is useful only if some practical pro- 
cedure for assembling such states can be devised. For the 
simplest possible quantum computer, comprising a single 
qubit, the process is trivial, as the thermal equilibrium 
density matrix has exactly the desired form, but with 
larger systems the situation is more complicated. For a 
system of n spin-^ nuclei, all of the same nuclear species 
(a homonuclear spin system), the 2™ eigenstates will be 
distributed across an evenly spaced ladder of n+ 1 groups 
of energy levels, with the number of (nearly degenerate) 
states within each group given by Pascal's triangle; the 
population of each state will be determined by the Boltz- 
mann equation. If the system contains several different 
nuclear species (a heteronuclear system), the situation is 
similar but slightly more complicated. Normally NMR 
experiments are conducted in the high temperature limit 
(see below), and so the pattern of excess populations will 
simply be proportional to the energy of each state. 

Assembling a pseudo-pure state from such a complex 
mixture might seem difficult, but it is in fact a fairly 
conventional problem from the viewpoint of NMR 0. 
The apparent problem is that reaching even a pseudo- 
pure state requires a non-unitary process, and the most 
obvious such process (relaxation to thermal equilibrium) 
leads to a state which is neither the desired state, nor 
unitarily related to it. Thus some additional non- unitary 
step is required. 

In fact there are two elements commonly used in NMR 
pulse sequences with non-unitary effects: magnetic field 
gradients and phase cycling. The first of these relies on 
the fact that the sample forms a macroscopic ensemble; 
by applying Hamiltonians which vary over the sample the 
ensemble averaged evolution can be non-unitary. This is 
most commonly achieved by momentarily destroying the 
spatial homogeneity of the main magnetic field (a Bq field 
gradient pulse |p5[ ), but similar effects can be achieved 
using spatially inhomogeneous RF fields (a B\ field gradi- 
ent pulse (3(| ) . The second approach relies on combining 
the results of several subtly different NMR experiments 
by post-processing; as this is done by classical methods, 
such processing is not confined to unitary transforma- 
tions. In conventional NMR this is referred to as phase 
cycling f3?| , and plays a central role in many experiments, 
although for many purposes it has now been replaced by 
the use of gradient techniques. 

Both of these approaches have been used to assemble 
pseudo-pure states. The original approach of Cory et al. 
0, based on magnetic field gradients, is in many ways 
the most satisfying, but an alternative "temporal averag- 
ing" scheme based on phase cycling Q has also proved 
extremely popular. Recently Knill et al. have described 



a simple and general approach |39| for building pseudo- 
pure states in a system of any size; their scheme may be 
used with either gradient or phase cycling techniques. 

In addition to these schemes, another quite different 
approach, logical labeling, was suggested very early on 
by Gershenfeld and Chuang [|| . Their approach is based 
on the observation that while the thermal equilibrium 
spin density matrix for an n spin system (n > 1) does 
not have the desired form, equation (|l|), subsets of en- 
ergy levels can be chosen which do have the correct pat- 
tern of populations; the computation is then performed 
within this subset of states. Although theoretically ele- 
gant this scheme appears to be more complex than the 
other approaches, and only two experimental implemen- 
tations have been reported Gol Ell . 



Scaling the system up 

While pseudo-pure states provide a practical approach 
for building small NMR quantum computers it is not pos- 
sible to simply scale this approach up to larger systems. 
The basic problem is the size of e, or rather the manner in 
which e scales with n, the number of spins in the system. 

The exact value of e will vary with details of the scheme 
used to prepare the pseudo-pure state, and so it is more 
useful to consider an upper bound, which gives the max- 
imum amount of pseudo-pure state which can possibly 
be extracted from thermal equilibrium. This is equal to 
the population difference between the lowest and highest 
energy states [^5| , which for an n spin homonuclear spin 
system is given by 



2smh(nhv/2kT) 
2™ cosh n (hv/2kT) 



(3) 



where hv is the Zeeman splitting. In the high tempera- 
ture limit (hv <C! kT) this expression simplifies to 



nhv/kT 



(4) 



Thus, in the high temperature limit the amount of 
pseudo-pure state which can be obtained decreases ex- 
ponentially with the size of the spin system. To over- 
come this it is necessary to use an exponentially large 
sample, or some equivalent approach such as repeating 
the experiment an exponentially large number of times. 
It therefore seems that the pseudo-pure state approach 
does not scale. 

In passing it should be noted that this is by no means a 
feature unique to NMR quantum computation: it occurs 
for any ensemble quantum computation scheme in the 
high temperature limit. Physically this is because such a 
system has 2™ levels, and the population deviations must 
be distributed among them; hence the excess population 
in any one state is inevitably exponentially small. 
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The exponentially small size of e has caused some au- 
thors to dismiss NMR quantum computation as a practi- 
cal approach; while this point of view has merit it may, as 
discussed below, be over hasty. More recently, it has been 
suggested that NMR might not be quantum mechan- 
ical at all! As NMR experiments are conducted in the 
high temperature limit, the density matrix is close to a 
maximally mixed state, and such high temperature states 
can always be decomposed as a mixture of product states 
(that is, states containing no entanglement between dif- 
ferent nuclei). As NMR states can be described without 
invoking entanglement, they can therefore be modeled 
classically, although the classical models involved may 
be somewhat contrived. While this conclusion is clearly 
correct it has proved difficult to develop classical models 
which fully describe NMR experiments 0] , and the real 
significance of these observations remains unclear. 

High fields and low temperatures 

As the problem with using pseudo-pure states arises 
from operating in the high temperature limit, the most 
obvious solution is to use either low temperatures or high 
fields, so that this limit no longer applies. Unfortunately 
NMR lies so far into the high temperature regime that 
this approach is unlikely to lead to success. 

The critical fields and temperatures required are given 
by ftrfB ~ kT; for 1 H nuclei and a magnetic field strength 
of 21.1 T (the largest NMR magnet currently available) 
this corresponds to T ~ 0.043 K. Reaching such temper- 
atures is possible, but clearly the sample will no longer 
remain in the liquid state! At such low temperatures 
only solid state NMR is possible. Alternatively if the 
sample is held at room temperature then a magnet with 
a field strength B ~ 150000 T will be required; this lies 
far beyond anything which is likely to be achieved in the 
foreseeable future. 

The arguments above do not entirely rule out the possi- 
bility that some combination of high fields and low tem- 
peratures might one day be used to achieve reasonable 
polarizations, and thus interesting pseudo-pure states (in 
particular, it may be possible to generate high polariza- 
tions under one set of conditions and then observe them 
in another), but in the short term this does not seem 
a particularly sensible approach. There are, however, a 
wide range of alternatives. 

Optical pumping 

A more subtle approach to increasing spin polarization 
is to use techniques such as optical pumping; this has the 
effect of decreasing the apparent temperature of the spin 
system without affecting the rest of the sample. Optical 
pumping techniques are in fact quite widely used within 



NMR, but it is not yet clear whether they can be usefully 
applied within NMR quantum computation. 

The best known optical pumping process within NMR 
is the spin-exchange optical pumping of noble gas nu- 
clei [Q, most notably the spin-^ nuclei 3 He and 129 Xe 
and the quadrupolar nucleus 131 Xe. The process involves 
conventional optical pumping of the electron spin states 
of alkali metal atoms, followed by spin exchange in bi- 
nary collision pairs (He) or van der Waals complexes 
(Xe). The resulting highly polarized noble gases have 
been used in a variety of NMR experiments, including 
NMR imaging Unfortunately noble gases are un- 

suitable for constructing NMR quantum computers as 
they consist of isolated atoms, that is systems containing 
only a single spin. It is in principle possible to transfer 
the polarization from the noble gas to more interesting 
species using a variety of cross polarization techniques 
p6| , [47| , f48[ , although it has so far proved difficult to 
obtain high transfer efficiency except when transferring 
polarization to surface species in microporous materials. 

The second common form of optical pumping in NMR 
is quite different: optical pumping in bulk semiconduc- 
tors such as Si, GaAs and InP |^] . This approach, which 
is related to the DNP schemes described below, is con- 
fined to solid state systems, and it is difficult to see how 
it might be used to improve current liquid state imple- 
mentations of NMR quantum computers. 

Other approaches 

There are many other techniques which can be used 
to increase the initial polarization in NMR experiments: 
the low sensitivity of NMR is perhaps its biggest draw- 
back in conventional spectroscopic studies, and seeking 
to improve sensitivity is a common research topic. 

Perhaps the most important technique for sensitivity 
enhancement is the nuclear Overhauser effect (NOE), 
which arises from the correlated relaxation of two or more 
nuclei. If the polarization of one nucleus is perturbed 
from its equilibrium value cross relaxation will transfer 
some of this perturbation to other nearby nuclei. This 
technique is widely used, both to enhance the polariza- 
tion of low sensitivity nuclei, and to probe internuclear 
distances [ p0| . The maximum polarization gain which 
can be achieved, however, is proportional to 75 /jj , where 
7s and 7/ are the gyromagnetic ratios of the sensitive and 
insensitive nuclei respectively, and so this method cannot 
be used to increase the polarization of 1 H, which has the 
highest gyromagnetic ratio among stable nuclei. 

A better approach for 1 H nuclei is to use the origi- 
nal Overhauser effect J5l], Q , which transfers polarization 
from electrons to nuclei. Initially greeted with scepticism 
p3[ , Overhauser's theoretical predictions were confirmed 
by Carver and Slichter [ f54| . who demonstrated huge 
Overhauser enhancements in the spectrum of metallic 
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7 Li. The Overhauser effect, and related phenomena col- 
lectively known as Dynamic Nuclear Polarization or DNP 
p5| , |56[ [57| , can be used to generate quite large polariza- 
tion enhancements in a range of solid systems contain- 
ing unpaired electrons, and when combined with optical 



techniques for pumping the electron polarization 49, 5q| 
dramatic enhancements can be observed p9[ . However 
the technique performs poorly in the liquid state. 

Another technique which gives enhanced polarizations 
is Chemically Induced Dynamic Nuclear Polarization, or 
CIDNP @. Despite the name CIDNP is unrelated 
to DNP; instead the non-equilibrium polarizations arise 
from a spin sorting mechanism which takes place dur- 
ing chemical reactions. This intriguing effect has proved 
a powerful tool for investigating a range of biomolccu- 
lar systems |30[ |Sl| , but the polarizations achievable are 
unlikely to be useful for quantum computation. 

Yet another possible approach is a family of exper- 
iments using para-hydrogen induced polarization, or 
PHIP jo^ ]. When hydrogen molecules are cooled into 
their rotational ground state, the Pauli principle dictates 
that the nuclear spin wavefunction of the two *H nuclei 
must be antisymmetric, and so the two nuclei must have 
opposite spin states. If para-hydrogen is used in an ad- 
dition reaction, for example adding H2 across a carbon- 
carbon double bond, then the product of the reaction 
will also have non-equilibrium spin states; these can be 
converted into a greatly enhanced polarization by con- 
ventional NMR pulse sequences. As para-hydrogen has a 
pure nuclear spin state it should in principle be possible 
to produce completely polarized molecules; in practice 
the enhancements are usually somewhat smaller [ |32| . 

Unlike many of the schemes discussed above, PHIP 
works well in the liquid state using fairly "normal" or- 
ganic molecules, and so may prove useful in NMR quan- 
tum computation. Unfortunately the scheme only allows 
the production of molecules which are highly polarized 
at two sites, and while it is in principle possible to use 
two or more addition reactions within the same molecule 
to produce polarization at four or more sites creating an 
entire spin system by PHIP addition reactions does not 
seem a plausible process. However, as discussed below, 
in some cases it may suffice to produce high polarizations 
at a single spin; in this case PHIP may indeed turn out 
to be a useful approach. 



Reinitialization 

One important point, frequently neglected in discus- 
sions of this kind, is that the methods described above 
are essentially methods for initializing an entire quantum 
computer at the start of a computation. They do not ob- 
viously permit the selective reinitialization of individual 
qubits in the middle of a computation. This makes it 
difficult to implement effective error correction schemes, 



as discussed in section below. 



Computational solutions 

In addition to the physical approaches outlined above, 
two computational approaches might allow the problem 
of low spin polarization to be bypassed. 

The first approach, due to Schulman and Vazirani 
p3| , uses computational methods to purify mixed states. 
Their scheme works by concentrating the polarization 
from a large number of weakly polarized spins into a small 
number of spins which become strongly polarized; alter- 
natively the calculation can be though of as occurring 
inside a small effectively pure subspace within the large 
spin space. Unfortunately the size of the subspace which 
may be extracted is 0(e 2 n), and so an n spin quantum 
computer requires a spin system with 0(n/e 2 ) spins; for 
the values of e achievable by direct cooling the result- 
ing overhead is enormous. If, however, the polarization 
can be substantially increased by other means then this 
might provide a good method for further purification. 

Another recent theoretical discovery may allow all 
these problems to be sidestepped, as it may be possible 
to perform many important quantum computations with- 
out starting from a pure (or even pseudo-pure) ground 
state. It has been known for some time that some quan- 
tum computations can be performed with a starting state 
comprising one pure qubit together with a number of 
qubits in the maximally mixed state ]64[ ] ; more recently it 
has been shown that Shor's quantum factoring algorithm 
can be performed in this way p5| (but see section be- 
low). If this proves to be a useful approach then it could 
simplify the problem of constructing an NMR quantum 
computer, as single qubit pure states may be substan- 
tially easier to reach than their multi qubit equivalents 
(for example, by using para-hydrogen based schemes). 



GATES 

In comparison with initialization, building quantum 
gates is a relatively simple process for NMR quantum 
computers. Indeed, as described above, many conven- 
tional NMR pulse sequences used in the molecular sci- 
ences can be viewed as sequences of quantum gates. 

Any quantum gate can be built out of one qubit and 
two qubit gates p6| , and early work concentrated on 
building these gates in two qubit NMR quantum comput- 
ers (NMR systems containing two spins). Unlike many 
other proposed implementations of quantum computa- 
tion, these gates cannot simply be transferred to larger 
spin systems without modification, and there was some 
initial concern that this transfer might prove difficult, 
involving an exponential increase in the complexity of 
the gate design. In fact, while the gates do need to be 
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modified, the transfer can be done fairly simply, and ap- 
proaches are known for doing this with only quadratic 
overhead. More recently some authors have moved on 
from simply imitating gates used in "conventional" quan- 
tum computers, and started designing gates which har- 
ness the full power of the NMR Hamiltonian directly. 



One qubit gates 

One qubit gates are easy to implement in NMR as 
they correspond to rotations of a single spin within its 
own subspace. This is most simply described using the 
Bloch vector model || ^ in which the state of a spin (or 
indeed any qubit) is represented by a point on the sur- 
face of the Bloch sphere. One qubit gates correspond to 
rotations on this sphere, and can be achieved using res- 
onant radiofrequency (RF) fields; the power and length 
of the RF pulse determines the rotation angle, while the 
phase of the RF radiation corresponds with the rotation 
axis within the xy plane. Rotations around other axes 
can be achieved by composite rotations, implemented by 
applying several RF pulses in sequence || . 

This description assumes that it is possible to apply RF 
pulses selectively to individual qubits. In more conven- 
tional proposals for implementing quantum computation 
such qubit selection is achieved using spatial localisation: 
each qubit is stored on some physical object with a well 
defined location. This approach is not possible within 
NMR for three reasons. Firstly, each qubit is not stored 
on a single spin, but rather in an ensemble of spins dis- 
tributed throughout the sample. Secondly, as the sample 
is in the liquid state individual spins are in continual 
rapid motion. Thirdly, the wavelength of RF radiation 
(around 1 metre) is huge compared with the separation 
between spins, rendering conventional spatial localisation 
impossible (although this could in principle be overcome 
by using techniques from magnetic resonance imaging 
(MRI) combined with enormous field gradients |38[ |69|). 

Instead of using spatial selection, NMR quantum com- 
putation relies on frequency selection to pick out individ- 
ual qubits. Different nuclear species have different gyro- 
magnetic ratios, and thus different resonance frequencies; 
therefore RF pulses on resonance with one spin will have 
little or no effect on other spins. With two or more nuclei 
of the same type the situation is slightly more complex, 
but even these will have different resonance frequencies 
as different chemical environments will result in differ- 
ent chemical shift (shielding) interactions fl. In this 
case it is necessary to use low power selective RF pulses 
("soft" pulses) or their multipulse equivalents, Dante 
sequences [70| , in order to excite one spin while leaving 
other spins with only slightly different resonance frequen- 
cies untouched (the excitation bandwidth of an RF pulse 
depends on its power) . This approach works well in small 
spin systems, but will be difficult in large systems as dis- 



cussed below. 



Two qubit gates 

Two qubit gates are slightly more complex, as they 
require some sort of interaction between pairs of spins. 
In liquid state NMR quantum computers this is provided 
by scalar spin-spin coupling ||. The NMR Hamiltonian 
describing two weakly coupled spins is 



H = ujjI z + ujsSz + ^Jis^IzS z 



(5) 



where, following NMR practice, the spin operators are 
described using product operators [72|, and energies 
are written as multiples of h. The one spin angular mo- 
mentum operators I z and S z are simply scaled versions of 
the Pauli matrices (the use of I and S to describe two spin 
systems is traditional j73|); u>i and uis are the resonance 
frequencies (Larmor frequencies) of the two spins, and 
J is is the scalar spin-spin coupling (J-coupling), usually 
measured in Hertz. Note that this simplified Hamilto- 
nian is only correct in the weak coupling limit, that is 
when \2irJ\ <C — uis\- 

The most obvious way to implement a two qubit gate, 
such as the controlled-NOT gate, is to use soft pulse tech- 
niques to selectively excite one spin only when its neigh- 
bouring spin is in one of its two eigenstates [Q; these 
two transitions have different energies as they are split 
by spin-spin coupling. This direct approach, picking out 
one transition from the multiplet corresponding to tran- 
sitions of a single spin, is in fact identical to the old se- 
lective population transfer experiment |p4f , and so this 
is the oldest method for building quantum logic gates. 
More recently this approach has been used to build two 
and three qubit gates directly [jl?], fll], (75| . 

This direct approach is, however, relatively unpopular 
festead^ost jmthOTsj7|Jl U [To| |l3| p| p[ |r| [l| g 



[H g fn], [78], [re], m UH, EJETSEiH 



prefer to use multipulse NMR techniques |2|, 13, |4|, |6| to 
sculpt the Hamiltonian, equation (|5|), into a more suitable 
form. This is analogous to the replacement of selective 
population transfer by the INEPT pulse sequence [ p2[ . 
The basic idea behind Hamiltonian sculpting is to use 
spin echoes || |9l], [)2j to refocus specific interactions in 
the Hamiltonian, thus creating an effective Hamiltonian 
obtained by rescaling elements of the original Hamilto- 
nian, or indeed completely deleting them. While this ap- 
proach to implementing two qubit gates can be described 
in a variety of different ways, they are all fundamentally 
equivalent. The central feature is the implementation of 
a controlled phase gate [^l) , such as 



4> 
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\fl e^J 



(6) 
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This can be decomposed in product operator notation as 

= cxp H x H x {-{hE) + I Z + S Z - 2I Z S Z )\ . (7) 

The first term {hE) can be ignored as it simply corre- 
sponds to an (undetectable) global phase; the remaining 
three terms can be achieved by sculpting the Hamilto- 
nian, equation (JsJ) , into the desired form. 

More recently a third approach for implementing two 
qubit gates has been suggested Q. Once again this 
method relies on the use of controlled phase gates, but 
the phase shifts are generated not by using conventional 
dynamic phases, but instead by geometric phases [p4| , 
such as Berry's phase Berry phases have been 

demonstrated in a wide variety of systems |94f| , includ- 
ing NMR Bq, j97f and the closely related technique of 
NQR |H II |0g, pal . They can be used to imple- 
ment controlled phase shift gates in NMR systems |)3| , 
but it seems that this approach has few advantages for 
NMR quantum computation over the more conventional 
dynamic approach; the idea may, however, prove useful 
in other systems [102 . 



Gate times 

It is not sufficient simply to show that a gate can in 
principle be implemented; it is also important to consider 
how long it takes to implement it. As discussed in section 
below, what matters is not the absolute time taken, but 
how this compares with the natural decoherence time. It 
is, therefore, important to consider what factors limit the 
rate at which NMR quantum, logic gates operate. 

One qubit gates are not only simple to build, they are 
also rapid to operate. As single qubit gates correspond 
to rotations, a reasonable measure is provided by the in- 
verse of the time required to perform a 2ir rotation. For 
a fully heteronuclear spin system (that is, a spin system 
containing only one spin of any given nuclear species) 
this rate is limited only by the available RF power and 
the breakdown voltages of the RF coils; typical values lie 
in the range 10-100 kHz. For a homonuclear spin system 
the situation is more complex as it is necessary to use se- 
lective excitation. In this case the gate rate is constrained 
by the frequency difference between the resonances which 
must be excited and those which must be left untouched. 
This can be seen either by considering a Fourier spectral 
model of excitation [Q, or by analysis of the "jump and 
return" pulse sequence [ 103 1 . Thus in homonuclear spin 
systems the one qubit gate rate is often below 1000 Hz. 

Two qubit gates are typically much slower. If they are 
implemented by direct selective excitation, then the sit- 
uation is the same as for one qubit gates in homonuclear 
systems, except that the relevant frequency splitting is 
the scalar coupling between the spins, Jis- A similar ar- 
gument applies if the gates are implemented using mul- 
tiple pulse techniques: the relevant rate is the inverse of 



the time required to achieve the "antiphase" condition, 
1/2 Jis- Scalar couplings are quite variable, but they are 
frequently less than 10 Hz; thus the gate times required 
for two qubit gates can be quite long. In our two qubit 
NMR quantum computer based on cytosine, a two qubit 
gate takes about 70 ms. If implemented using geometric 
phases 93 , the time required is even greater. 

In order to reduce these two qubit gate times it is nec- 
essary to increase the size of the spin-spin coupling con- 
stant. In general this is not possible, as J/s is fixed by 
the chemical system chosen (it is, of course, possible to 
choose a new system with a larger coupling constant, but 
this approach is obviously quite limited). It is, however, 
possible to change the apparent size of J/s by partially 
reintroducing dipolar couplings. The scalar coupling is in 
fact a quite small effect in comparison with the principal 
spin-spin coupling interaction, dipolar coupling. This 
arises from the direct magnetic interaction between two 
magnetic dipoles, and in the high field approximation jl]] 



Mo7/7sfr(! ~ 3 cos 2 6) 
8nr 3 



(3I Z S Z -I-S), (8) 



where r is the length of the internuclear vector r, and 9 
is the angle between r and the magnetic field. In solid 
samples this results in coupling between one spin and 
all its neighbours, but in liquids and solutions rapid mo- 
tion causes r, and thus 8, to fluctuate, so that Hd is 
replaced by its isotropic average, which is zero. For this 
reason dipolar coupling has little direct effect in liquid 
state spectra (it cannot be entirely neglected, as it is one 
of the main sources of spin relaxation) . The scalar cou- 
pling is a correction to the dipolar coupling which arises 
from the Fermi contact interaction: valence electrons can 
interact with two or more nuclei and thus mediate an 
interaction between them. Like dipolar coupling scalar 
coupling is anisotropic, but its isotropic average is non- 
zero; thus the relatively small isotropic scalar coupling is 
the dominant spin-spin interaction in the liquid state. 

Dipolar coupling is removed because isotropic tum- 
bling reduces it to its isotropic average. If the motion is 
anisotropic then the average may become non-zero; this 
can occur if the anisotropy of the molecular magnetic 
susceptibility causes molecules to align slightly with the 
magnetic field, resulting in small residual dipolar cou- 
plings. The effect is largest when the magnetic suscepti- 
bility is high, and this approach has proved helpful in 
investigating proteins bound to oligonucleotides [ 104 1 . 
Larger alignments, and thus larger residual couplings, 
can be observed using liquid crystalline solvents or cosol- 
vents [ 105 , which align strongly in the magnetic field and 
then themselves act to order other dissolved molecules. 



This approach has been used 1 106 to implement 



Grover's quantum search on a two qubit NMR quantum 
computer based on chloroform dissolved in a liquid crys- 
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tal solvent. The effective Hamiltonian in this system is 
H = + u s S z + tt( J is + 2D)2I Z S Z (9) 



where D is the residual dipolar coupling [106], allowing 
the gate rate to be increased by a factor of 8. It is not 
yet clear how generally useful this approach will be. 



Two qubit gates in larger systems 

The description above is not only confined to one and 
two qubit gates; it is in fact confined to such gates in 
two qubit computers. If these gates are to be used to 
build networks with three or more qubits it is necessary 
to consider how they will function in larger spin systems. 

The situation for one qubit gates is relatively straight- 
forward, as these can still be implemented using selective 
RF pulses (although the problem of frequency selection, 
discussed below, becomes more serious). For two qubit 
gates, however, a more detailed analysis is necessary. For 
a general n spin system, the Hamiltonian is 



H 



i<k 



nJik2I l I 



(10) 



including n Larmor frequency terms for each spin, and 
a total of n(n — l)/2 spin-spin coupling terms, connect- 
ing every pair of spins. In real systems some of the Jik 
coupling constants will be so small that they can be ne- 
glected, but it is still useful to consider the most general 
case. Consider, for example, a system of three spins, 
conventionally called I, R and S; in this case there will 
be three Larmor frequency terms, and three spin-spin 
couplings, J IR , Jis, and J RS . 

It is in principle relatively simple to implement gates 
using the direct method. Each transition is split into four 
(under the influence of two couplings); if selective exci- 
tation is applied at just one of these four frequencies this 
corresponds to a doubly controlled three qubit gate, such 
as the Toffoli gate. Alternatively by exciting two tran- 
sitions two qubit gates, such as controlled-NOT, can be 
achieved. In practice in large spin systems it will become 
difficult to pick out the right set of transition frequencies, 
and it is likely that some transitions will overlap. 

The more common approach is to use Hamiltonian 
sculpting to convert equation ( |To| ) into the desired form. 
Clearly this requires refocusing the n — 2 additional Lar- 
mor frequency terms, as well as the (n — l)(n — 2)/2 
extra spin-spin couplings. This can be achieved using 
spin echoes J82], p8| . Unfortunately it is not possible to 
use just a single echo to refocus all the interactions; in- 
stead it is necessary to consider how the echo sequences 
interact with one another. The simplest method is to 
nest spin echoes applied to each spin within one another, 
but this naive approach requires an exponentially large 
number of refocusing pulses. Fortunately this problem 



can be si destepped by using efficient refocusing sequences 
p8| , 107 1, which allow refocusing to be achieved with a 
quadratic overhead. 

It is rare to find a large spin system where all the cou- 
plings have significant size. Coupling is a fairly local 
effect and so large coupling constants are only found be- 
tween spins that are fairly close within the spin system; 
thus the coupling network is described by a non-complete 
graph |Q . This greatly simplifies the problem: not only 
does it reduce the number of couplings which have to be 
refocused, but it also simplifies the echo patterns needed 
to refocus the Larmor frequencies fl 



Thus the over- 
head required for refocusing echoes is greatly reduced. 

It might seem that it would not be possible to imple- 
ment all the desired logic gates in such a spin system, as 
some of the requisite spin-spin couplings are missing. In 
fact this is not a problem, as long as every pair of spins 
is connected by some chain of couplings. Quantum SWAP 
gates |8C], |83j can be used to move quantum information 
along this chain; the overhead imposed by these SWAP 
gates is at worst linear in the size of the spin system. In 
most cases the advantages of using a partially coupled 
spin system significantly outweigh the disadvantages. 



Multiqubit gates 

As described above, NMR techniques can be used to 
implement conventional one and two qubit gates, and 
thus any desired quantum network. This approach, im- 
itating existing theoretical models of quantum comput- 
ers, was adopted by all the early papers on NMR quan- 
tum computation. It is not clear, however, that this is 
the best approach: it might be more sensible to consider 
what types of gates NMR systems are good at providing, 
and then seeing whether these are valuable in quantum 
computation. To take a simple example the natural two 
qubit gate for NMR systems is the controlled phase shift, 
not the controlled-NOT gate, and in many quantum al- 
gorithms phase shifts are exactly what is needed. For 
example, NMR quantum computers can readily imple- 
ment Grover's quantum search without using an ancilla 
qubit to convert oracle calls into phase shifts |lj, |l5|, [l6) . 

A more complex, but more interesting, example is the 
implementation of multiqubit gates. As the NMR Hamil- 
tonian, equation (10), contains terms connecting multiple 
pairs of qubits, it should be possible to use this Hamil- 
tonian to build certain multiqubit gates directly. This 
was briefly discussed above, when considering the direct 
implementation of Toffoli gates in a multi-spin system, 
but Hamiltonian sculpting should allow the approach to 
be used more widely. To date this idea has received only 
brief analysis p7| , but this is likely to be a productive 
source of gates in the next few years. 
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The problem of selective excitation 



DECOHERENCE 



From the descriptions above it might seem that im- 
plementing quantum logic gates in NMR quantum com- 
puters is essentially solved, and that the solutions scale 
fairly well. This is almost, but not quite, true as one 
major problem remains: selective excitation. 

All the techniques described above rest on the assump- 
tion that it is possible to address individual qubits, so 
that interactions can be refocused. For NMR quantum 
computation (or indeed any NMR pulse sequence) this 
is achieved using frequency selection, rather than more 
conventional spatial localisation techniques, and this ap- 
proach is only possible if the NMR transitions all have 
well separated frequencies; in particular it is simplest if 
the separation between any pair of Larmor frequencies 
is much greater than the width of the NMR multiplets, 
that is the sum of the spin-spin couplings. Note that 
in contrast with conventional NMR experiments it is not 
sufficient simply to selectively excite one spin; it is also 
essential that other spins remain completely unaffected. 

With small spin systems this is fairly easy to achieve, 
but with larger systems it can become quite difficult. For 
example, the range of Larmor frequencies found for *H 
nuclei in simple organic compounds, and working at a 1 H 
frequency of about 500 MHz is only about 5000 Hz, and 
this limited frequency range can soon "fill up" . Increas- 
ing the number of spins not only increases the number 
of NMR multiplets, but also increases the width of each 
multiplet by introducing more spin-spin couplings. With 
other nuclei the situation is similar, although usually less 
serious (the X H frequency range is unusually narrow). 

A partial solution is provided by turning to heteronu- 
clear spin systems. As the NMR frequencies of differ- 
ent nuclei are very different, it is trivial to achieve nu- 
cleus selective, and thus spin selective, excitation. The 
comparative simplicity of heteronuclear NMR quantum 
computation explains its early, and enduring, popular- 
ity. Unfortunately this approach can not be continued 
indefinitely, as the number of suitable nuclei is small, the 
obvious candidates being 1 H, 13 C, 15 N, 19 F and 31 P. 

It seems likely that the problem of selective excitation 
will prove a serious difficulty in constructing large NMR 
quantum computers. Although it is difficult to assess ex- 
actly what the limit will be, it is notable that the largest 
number of spins of one nuclear type used to date is six X H 
nuclei [ j35| ; one other paper has described computations 
involving four 13 C nuclei and three X H nuclei |3{|, while 
all other authors have used at most three spins of any one 
nuclear type. Assuming that it is practical to address six 
spins of each of the five nuclei listed above, this suggests 
a limit of around 30 qubits imposed by the problem of 
selective excitation. Actually designing and synthesising 
such a spin system is another problem entirely. 



In order to perform large quantum computations it 
is essential that errors be kept under control. In prac- 
tice this means that the decoherence time must be very 
long in comparison with the gate time, although meth- 
ods of error correction (see section ) allow this criterion 
to be slightly relaxed. The situation for NMR quantum 
computation might appear extremely good, as systems 
with very long relaxation times are known: for exam- 
ple, the spin-lattice relaxation time (Ti) of 129 Xe can be 
thousands of seconds ]108| . As scalar spin-spin coupling- 
constants can reach hundreds of Hertz, and dipolar cou- 
plings can be even larger, a naive calculation suggests 
that it should be possible to implement about 10 4 -10 6 
gates before decoherence becomes a serious problem. 

In fact this calculation is meaningless for two reasons. 
First, these extremely long relaxation times are always 
spin-lattice relaxation times; the spin-spin relaxation 
times (T2), which provide a better measure of the deco- 
herence rate, are much shorter (typically below ten sec- 
onds). Secondly, the relaxation times and gate rates are 
taken from different spin systems, and it is not possi- 
ble to simultaneously achieve them in a single molecule. 
The long relaxation times observed for 129 Xe gas arise 
precisely because Xe gas atoms have only weak interac- 
tions with their neighbours, and such systems are of little 
apparent use for quantum computation as they provide 
no mechanism for logic gates. As the scalar coupling is 
related to an underlying dipolar coupling, and dipolar 
coupling is one of the two principal sources of relaxation 
in spin-i nuclei, spin systems suitable for constructing 
NMR quantum computers inevitably possess shorter re- 
laxation times. 

One further feature of decoherence in NMR quantum 
computation deserves consideration. As NMR systems 
are ensemble quantum computers, the effect of decoher- 
ence is not simply to introduce errors, as occurs with 
more conventional designs. Rather these errors must be 
averaged over the ensemble, and if, as is often the case, 
these errors are fairly random in character, the overall ef- 
fect is that the error terms will largely cancel out. Thus 
the principal effect of decoherence is to reduce the ap- 
parent signal strength. This is clearly visible in some 
quantum counting experiments [119, 109 , where decoher- 



ence appears as an exponential decay in signal. 

It is difficult to estimate a realistic limit to NMR 
quantum computation arising from decoherence; how- 
ever, current experiments have been performed involving 
hundreds of logic gates [jl9L |109| , and it seems likely that 
the limit is about one thousand gates. 
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QUANTUM ERROR CORRECTION 

Quantum error correction, and its companion fault tol- 
erant computation 11C, 111, 112, 113, 114, 115], play a 
central role in considerations of the practicality of large 
quantum computers: while small quantum computers 
may handle errors fairly well the fragility of highly entan- 
gled states renders large quantum computations highly 
vulnerable. Error correction tackles this by diagnosing 
these errors and correcting them, while fault tolerant 
computation provides methods for minimizing the spread 
of errors, and in particular reducing the impact of errors 
in error correction schemes. Simple examples of such 
schemes have been implemented on NMR quantum com- 
puters 0, |§. 

The importance of these techniques can hardly be over- 
stated; until their discovery many authors believed that it 
would be completely impossible to build a large quantum 
computer as it would be hopelessly error prone. Unfor- 
tunately this benefit comes at a price: error correction 
involves a substantial overhead as each logical qubit must 
be encoded using many ancilla qubits. While thi s ov er- 
head has been substantially reduced by new codes 115] it 
is still at least an order of magnitude. More importantly, 
however, this assumes that it is possible to reuse ancilla 
qubits in order to repeatedly correct errors; in turn this 
requires the ability to reinitialize ancilla qubits at will. 
With current NMR quantum computers this cannot be 
achieved, and it is necessary to use a supply of fresh 
ancilla qubits at each stage. In this case the overhead 
becomes so large as to be completely impractical. 



READ-OUT 

Once a quantum computation has been performed it 
is necessary to use some read-out scheme in order to ex- 
tract the result. As NMR quantum computation is imple- 
mented using an ensemble of spin systems, read-out in- 
volves ensemble measurements, and thus expectation val- 
ues. This can have a number of profound consequences. 

In simple cases NMR read-out is little different from 
more conventional schemes. Suppose that a quantum 
computation ends with the answer qubits in eigenstates. 
In this case it is only necessary to determine whether 
each qubit is in state |0) or |1), which is equivalent to 
determining the expectation value of a z . This can be 
achieved either by exciting the corresponding spin with 
an RF pulse and observing the phase of the NMR signal 
JTo| , [l9| , or by examining the multiplet structure in the 
NMR spectrum of a neighbouring spin, or by a combina- 
tion of these techniques |l3| . 

A more interesting situation occurs when the algorithm 
ends with one or more answer qubits in superposition 
states. A conventional quantum computer will return 
one of the corresponding answers at random, while an 



ensemble quantum computer will return an ensemble av- 
erage over the set of all possible answers . In most 
cases this result is not particularly useful, and it is neces- 
sary to recast the algorithm so that a single well defined 
result is obtained (|, In some special cases, however, 
ensemble measurements can be advantageous [ pTsj ] - 

A similar situation arises when NMR techniques are 
used to implement phenomena such as quantum telepor- 
tation |7^, 116 1. Traditional teleportation schemes use 
strong measurements to project an unknown quantum 
state into the Bell basis; classical results from these mea- 
surements are then used to determine which of a set of 
unitary transformations must be used to finish the pro- 
tocol. In NMR teleportation, however, such projective 
measurements and classical readout cannot be used. In- 
stead it is necessary to use conditional evolution to per- 
form the final unitary transformation J79| . 

This lack of projective measurements may also have 
consequences for quantum computation using a single 
pure qubit |54| ^35| . The model of Knill and Laflamme |6^| 
assumes only ensemble measurements, and so accurately 
reflects the nature of NMR quantum computation, while 
the work of Parker and Plenio ] |65| assumes that projec- 
tive measurements can be made. It remains to be seen 
how significant this difference actually is. 



CONCLUSIONS 

It is useful to draw some overall conclusions about the 
potential usefulness of liquid state NMR as a technique 
for implementing large quantum computers. For small 
quantum computers liquid state NMR techniques are well 
ahead of the competition: indeed in most areas of ex- 
perimental quantum information processing there quite 
simply is no competition! However, there are several se- 
rious difficulties with extending this approach to large 
systems, and it seems unlikely that any very large liquid 
state NMR quantum computer will ever be built. 

The issue most commonly raised is initialization: the 
pseudo-pure states used in NMR are far from pure, with 
typical polarizations below 10 -4 , and, more seriously, an 
exponential fall off as the size of the computer is in- 
creased. This alone would appear to limit liquid state 
NMR quantum computers to about 30 qubits. 

In fact this assessment is probably too pessimistic. 
There are many techniques for increasing signal 
strengths, and while none of them offer an immediate 
solution, several have potential. Furthermore, recent the- 
oretical results suggest that pure states may be less im- 
portant than previously believed. If low polarizations 
were the only difficulty preventing the construction of 
a large scale NMR quantum computer, then it seems 
highly likely that a solution would be found. Unfortu- 
nately there are other, more serious, problems. 

Constructing quantum logic gates in small systems is 



11 



easy, indeed almost trivial; this simplicity explains the 
rapid initial progress in NMR quantum computation. 
Early concerns about a potential exponential increase in 
the complexity of implementing these gates have proved 
unfounded, with the discovery of methods for implement- 
ing two qubit gates in multi-spin systems with (at worst) 
quadratic overhead. Little thought has, however, been 
paid to the problem of selective excitation and the crowd- 
ing of frequency space; in my personal opinion this is 
likely to be the first serious barrier to building NMR 
quantum computers with many more than ten qubits. 

The problem of decoherence is, of course, common 
to all potential implementations of quantum computers. 
The relatively long decoherence times of nuclear spins, 
which appear to make them good candidates as qubits, 
arise because they interact only weakly with their envi- 
ronment. Unfortunately these weak interactions are also 
manifested as low gate rates, and so the ratio of decoher- 
ence time to gate time is not as large as one might hope. 
Nevertheless NMR decoherence processes are fairly well 
behaved, and are not likely to prove an insuperable prob- 
lem in systems of ten qubits. 

The difficulty of performing practical quantum error 
correction schemes, arising from the lack of selective 
reinitialization, is a serious problem for building large 
quantum computers; this, however, is only likely to be- 
come an issue if all other difficulties are solved. 

The inability to perform projective measurements, and 
the difficulty of dealing with ensemble averaged data, is 
probably not a serious problem in its own right. This 
property makes it difficult to use NMR experiments to 
test fundamental questions in quantum mechanics, but 
has few major implications for computation. Indeed in 
some cases, such as when considering the effects of de- 
coherence, it can actually be an advantage. Unlike some 
other issues, however, this limitations appears to be in- 
herent in the liquid state NMR experiment, and it is not 
at all clear how it can be bypassed. 

Finally, it is interesting to compare current liquid state 
NMR quantum computers with Kane's radically differ- 
ent, but ultimately related, proposal for a solid state 
NMR quantum computer [Q. It is notable that Kane's 
proposal keeps many of the advantages of NMR, but also 
manages to tackle some of the most serious difficulties 
described above. Firstly, Kane's proposal uses external 
"gates" to modulate both the Larmor frequencies and 
the spin-spin coupling constants of nuclei. This should 
remove all the difficulties with implementing quantum 
logic gates. Secondly, the proposal includes a scheme 
for making projective measurements on single spins, thus 
also providing a simple initialization scheme. 

Even if, as seems likely, a large liquid state NMR quan- 
tum computer is never built, quantum computation will 
still owe much to NMR. By providing the first work- 
ing quantum computers, no matter how small, NMR 
has reinvigorated the field. Tricks long known to NMR 



spectroscopists are now being applied in NMR quantum 
computations, and many of these will have applications 
in other technologies. NMR quantum computation has 
had, and still has, much to offer. 
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